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Abstract 

The explanation of the large CP asymmetries in — )■ [7r~^7r~ , K~) decays observed by the 
LHCb collaboration is likely to call for new physics beyond the CKM paradigm. We explore new 

■ contributions caused by the color-sextet scalar diquark, and demonstrate that the diquark with 

■ the mass of order 1 TeV and nominal couplings with quarks can generate the CP asymmetries 
r~| . at the percent level. Using the experimental data on branching ratios and CP asymmetries of 

—7- (vr+vr", K^K~), we derive the constraints on the diquark mass and couplings, which can be 
further examined in hadron colliders in the dijet final states. 
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It is generally anticipated that both direct and indirect CP asymmetries (CPAs) in the 
charm sector are quite small in the standard model (SM). Any observation of the large CPA 
in decays will presumably imply that the underlying theory is out of the scope of the 
SM. 

Recently based on the 0.62 fb~^ of data collected in 2011, the LHCb collaboration [l| 
has measured the difference between the time-integrated CP asymmetries in the decays 
L)0 ^ K+K- and tt+tt", AAcp = Acp{D^ K+R-) - Acp{D^ tt+tt-), given 

by 

AAcp = (-0.82 ± 0.21(stat.) ± 0.11(sys.))% , (1) 

where the first uncertainty is statistical and the second is systematic. The quantity 
Acp{D^ f) is defined as 

Acp{D ^ /) = rpo^j) + r(DO^/) ' 



with / = K^K , vr+vr . By contrast, results released by the CDF collaboration j2| based 
on 5.9 fb"^ of the integrated luminosity are somewhat less conclusive, given by 



Acp{D^ TT+TT") = (+0.22 ± 0.24 ± 0.11) 



Acp{D^ = (-0.24 ± 0.22 ± 0.09)% , (3) 

while the previous world averages from Heavy Flavor Averaging Group jsj] in 2010 are 

Acp{D° ^ TT+TT-) = (+0.22 ± 0.37)% , 
Acp{D° K+K^) = (+0.16 ± 0.23)% . (4) 
The new world average for AAcp from Eqs. ([T]), ([3]) and (jl]) is found to be [4| 

AAcp = -(0.645 + 0.180)%, (5) 

which is about 3.6cr away from zero. 

Contributions to AcpiD'^ — > /) contain both direct {A'^,p{D^ — )■ /)) and indirect 
{A^^p{D^ — )■ /)) parts, and from the LHCb report jl|] one has 

AAcP ~ AA'^jp + (9.8 + 0.3)%A^^ , (6) 

where AAfi'p = Ag^(L)° K+R-) - A'^J'p{D° tt+tt"), and A^^ = A^^{D° /), 
which is universal for / = K^K" and tt+tt" and less than 0.3% due to the mixing param- 
eters. Clearly, the LHCb data in Eq. ([1]) is dominated by the difference of the direct CP 
asymmetries, AA'^'^p. 



In order to have a nonzero direct CPA, two amplitudes Ai and A2 with both nontrivial 
weak phase difference 9w and strong phase difference Ss are called for, leading to 
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- 2r f sin 9w sin ^5 , (7) 



with Vf = \A2\/\Ai\. In the last step, a hierarchy of rj- <C 1 has been adopted. The SM 
description of the direct CPA for — / arises from the interference between tree and 
penguin contributions, in which decay amplitudes take the generic expressions 

= V:^Vu,TsM + KtKftPsM, (8) 

with q = s for / = K^K~ and q = d ioi tt^tt". Besides the hierarchy in the CKM matrix 
elements V*gVuq ^ Kfe^^f" penguin amplitudes are also suppressed by loop factors. Even 
in the limit Psm ~ Tsm, the ratio of the decay amplitudes is still very small rj ~ 0.0007, 
leading to a tiny CPA which is far below the central value in Eq. (|5]). As a result, if we took 
the data by the LHCb seriously, a solution to the large AAcp would be to introduce some 
new CP violating mechanism beyond the CKM. 

By neglecting the small SM penguin contributions, the decay amplitude of D° —>/(/ = 
K^K~ or TT+TT^) with new physics contributions could be parametrized as 

Aj = [T'sm (1 + Pe^^-) + E'sMe''^] , (9) 

where 9w is the new weak CP phase ranging from to tt, while p is associated with the 
new physics effect with an arbitrary sign, i.e., sign(p) = ±1. In the above equation, T'gj^ 
corresponds to the W emission diagram, while -E^^/ is from the annihilation type of the 
W-exchange diagram. We note that since the final state interactions make dominant con- 
tributions to the W-exchange diagram, without losing generality, we regard that the short- 
distance (SD) effect of new physics on such topology could be ignored. Consequently, in this 
circumstance new physics plays an important role for the emission topology. From Eqs. ([7]) 
and (ED, we find 



Adir iT^Q . n _ ^T'smE'smP sill sin 6s , . 



with a{() = T'gj^jC + Egjy.je^^^ and ^ = 1 + pexp{i9w)- If all quantities in the SM are under 
control, p and 9w are the only free parameters. 



Recently, a number of theoretical studies |5l-ll7j have been performed to understand the 
LHCb and CDF data. In this brief report, we would like to use the scalar diquarks as the 
sources of new physics. Although the introduction of scalar diquarks in the literature has its 



physical reason, such as solving the strong CP problem |18H20|. here we only explore their 
consequences in D-meson processes. The combinations of two quarks give many possible 
types of diquarks, such as (3,1,-1/3), (6,1,1/3), (3,3,-1/3), (3,1,-4/3), and (6,1,4/3) 
under the standard group of SU{3)c x SU{2)l x f/(l)y, among which the (3, 1, —1/3) and 
(6, 1, 1/3) diquarks are very interesting for avoiding the strong correlation in flavor couplings 



and the strict constraint from the tree induced D-D mixing [21|. In what follows we will 
concentrate on the color sextet = (6, 1, 1/3) to illustrate its effect on the direct CPAs in 
D^ —7- K~^K~, TT+TT", but a similar study could be applied to the color triplet boson. 
We first write the interaction between quarks and Hq as 



^He - fijdiaCPLU'jf^HQ^ 



h.c. 



(11) 



where fij denotes the couplings between the diquark and quark flavors, C = Z7 7 is the 
charge-conjugation matrix, Pl[r) = (1 T 75)/2 is the chiral projection operator, and Hq'^ 
is a weak gauge-singlet and colored sextet scalar with a and /3 being the color indices. In 
terms of the interactions in Eq. fllip . flavor diagrams for D decays are given in Fig. [TJ After 
integrating out the highly virtual diquarks, the corresponding interactions for c — > udd{ss) 
are derived as 
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with {q"q')v+A = q'Yil 
form factors, defined by 



01 = {uq)v+A{qc)v+A , 

01 = iUaqp)v+Aiqi3Ca)v+A , (12) 

75)g' and q = d,s. Based on the decay constants and transition 



{0\qYl5q\P{p))=^fpP^, 
{P{p2)\q7,c\D{p,)) = fr{e){Q, - y^k,} 
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(13) 



respectively, with Q = pi + P2 and k = pi — p2, the diquark contribution to D — )• / is given 
by 



pe 
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(14) 



where oi = cAu) + C2{p){l/Nc + and a[ 



coefficients 



{1 / Nc + xifJ')) are the effective Wilson 



22| with xil^) being the nonfactorizable contribution. In the large Nc limit, as 
the nonfactorizable effect could be simplified as x = — [23j, the nonfactorizable part 
will be smeared by the operator 



C/3 
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FIG. 1. Diquark-mediated flavor diagram for — t- (vr'^vr ,K~^K 



Before the numerical study, we discuss how to escape the stringent constraint from the 
D — D mixing. Using the results in Ref. 2l|, the D mixing parameter induced by box 
diagrams could be formulated by 



xd 



D 



15367r2 



rriH 



(15) 



Since the flavors (denoted by i) in the internal loops include d, s and b quarks, the constraint 
from xo could be released if a cancellation occurs in Yli ftcL^ 



In our analysis, the input values of the SM are taken as 
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T'sj^,j{TTn,KK) = (3.01,4.0) X 10"^ GeV, 

E'sm{t^t^.KK) = (1.3,1.6) X 10"^ GeV, 

^^(TrTr, KK) = (145, 108)° , ai = 1.21, 

Vus = -Vcd = 0.2252 , Vcs = 0.97345, = 0.97428 , 

m^^K) = 0.139(0.497) GeV , mo = 1-863 GeV, 



(16) 



where the resulting branching ratios (BRs) for D° — )> (vr vr"*", K K^) are estimated as 
(1.38,3.96) X 10-^ while the current data are B{D^ 7r-7r+) = (1.400 ± 0.026) x lO'^ 
and B{D^ K-K+) = (3.96 ± 0.08) x 10-^ {m]. Since f*Jdu and f^Jsu are independent 
free parameters, to simplify our calculation, we adopt two benchmark schemes: (1) f^^fdu ~ 
fscfsu and (II) fdcfdu ~ —fscfsu- Consequently, the involved parameters in the analysis are 
Ow = arg{f;jsu) and \f;jsu\/mjj. 

An estimate of the scalar diquark contribution is given as follows. The current measure- 
ment of the dijet cross section from the hadron collider puts the limit of the scalar diquark 



mass, see Ref. 
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26| for instance, 

niH, > 1.9TeV, (17) 

where a normal diquark-quark coupling is used. When the diquark decay is taken into 
account, this value may get reduced. Assuming the mass of order 1 TeV and normal couplings 
for the diquark, we find from Eq. f|T4|) that the NP contribution is at the percent level 
compared to the SM contribution. As a result, such a diquark is able to explain the large 
CPA data by the LHCb, while its effects to branching ratios will be up to a few percent. 

Since the involved parameters in Cabibbo allowed processes, e.g., D — > txK, are different 
from the singly Cabibbo suppressed decays, with the assumption of ^ ■ /*,/m — ^ 0, the BRs 
for — (tt+tt", fC+A'") are the potential constraint. Thus, we will take the data of the 
BRs with errors to constrain the free parameters. For Scheme I, we show the CPA difference 
defined in Eq. ([6]) as a function of = ±|/*^/su|/m^ and 6w in Fig. [2l The solid lines 
correspond to the upper and lower bounds of the LHCb results. The dashed (blue) and 
dash-dotted (red) lines denote the experimental BRs for the decays — {ti^ti~ ,K^K~), 
respectively. From these curves, we obtain the allowed region of the parameters as 

-2 X 10^^ <^ < -0.6 X 10"^ 

0.3 <ew < 2.25. (18) 

Similarly, the results in Scheme II are presented in Fig. |3l In this scheme, we find that 
the magnitude of AA^p cannot fit the data within la errors. The reason is that the direct 
CPA of — 7- 7r"''7r~ in scheme II has the same sign with the one of — K^K~ . In order 
to understand the effects on each direct CPA A^*p defined in Eq. ([7]), by using Scheme I, 
we display A^*p for — tt+tt" and K^K~ decays in Figs. and |4}d, respectively. The 
results, even the sign, are consistent with the current CDF data shown in Eq. ([3]). 
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FIG. 2. AA'^p (solid) in units of 10~^ as a function of ^ = ib|/*^/sM|/m|^ and 6^ in the diquark 
model, where the dashed (blue) and dash-dotted lines denote the data with errors of BRs for the 
decays — )• {7:~^ir~ , K~) in units of 10~^, respectively, while the shadowed region (green) is 
the combined constraint from BRs and direct CPAs. 
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FIG. 3. Legend is the same as that in Fig. [2l but the magnitude of AA'^,^ is smaller than the LHCb 
results of la. 

In summary, it is clearly an exciting moment if any CPA in D-meson decays is observed 
at the percent level as the SM prediction is far below 1%. Motivated by the recent LHCb 
measurement on the time integrated CPA in — )■ {tt'^tt' , K^K~) decays, which appears 
to be inconsistent with the SM result, we have studied the contributions of the colored 
scalar boson and used the color sextet (6, 1, 1/3) as the illustrator. In this diqaurk model, 
the serious constraint from the D mixing parameter induced by box diagrams could be 
avoided if a cancellation among different flavor couplings occurs. We have found that the 
the induced direct CPAs of — )■ {'n'^n~ , K~^K^) decays can fit the recent LHCb results 





and are consistent with the CDF current measurement. 
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